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Oh' Abstract. We discuss some of the main open issues in the evolution of Active Galactic 

Nuclei which can be solved by the sensitive, wide area surveys to be performed by the 
proposed Wide Field X-ray Telescope mission. 
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1 . Introduction across a galaxy's lifetime. Support for this hy- 

— u. ■ pothesis comes from several lines of evidence 

(SI ' The observed scalin g relations between the such as i} the match bet ween the mass func- 

O ■ structural properties of massive galaxies (bulge tion of BHs grown thr ough A GN phases and 

mass, luminosity, and stellar velocity disper- that observed in local gala xies (iMarconi et all 

sion) and the mass of the black holes (BHs) [20041 IShankarlTall 12004): ii) the cosmologi- 

at their centers suggest that galaxy assem- cal "downsizin g" of both n uclear activity and 

bly and black hole growth are closely re- starformation dUedaet al.ll2003l: ICowie et~all 

lated phenomena. This is often referred to [199^); iii) the feedback produced by the AGN 

as BH/galaxy co-evolution. The BH growth on to the galaxy interstell ar medium through gi - 

is thought to happen primarily through effi- ant mo i eC ular outflows (iFerudio et alJ kOlO). 

cient accretion phases accompanied by the re- A numb er of semi-analytic models (SAMs) 

lease of kinetic and radiative energy, part of have been p ropose d over about the past 

which can be deposited into the galaxy inter- decade to explain the BH/galaxv co-evolution 

stellar medium. Active Galactic Nuclei (AGN) dKauffmann & Haehneld 120001: iMonaco et alJ 

are then believed to represent a key phase |2007t iMenci et alJ 120081: iMarulli et alJ I2008F 

lLamastra et al.ll201(J) . These models follow the 
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across cosmic time, either through the Press- 
Schechter formalism or through N-body sim- 
ulations, and use analytic recipes to treat the 
baryon physics within the dark matter halos. 
A common assumption of these models is that 
mergers between gas-rich galaxies trigger nu- 
clear activity and star formation. Recently, a 
BH/galaxy evolutionary sequence associated 
to "wet" galaxy me rgers has been proposed 
dHopkins et al.l2008l) . in which an initial phase 
of vigorous star formation and obscured, possi- 
bly Eddington limited, accretion is followed by 
a phase in which the nucleus first gets rid of the 
obscuring gas shining as an unobscured QSO, 
then quenches star formation, and eventually 
fades, leaving a passively evolving galaxy. 

While being successful in many respects, 
this picture is still very partial and many fun- 
damental pieces are missing to get a satisfac- 
tory understanding of BH/galaxy co-evolution, 
such as the very first stages of this joint evo- 
lution (e.g. at redshifts z > 6), the cosmo- 
logical evolution of nuclear obscuration, the 
triggering mechanisms and environmental ef- 
fects of nuclear activity. In this contribution it 
will be shown how the proposed Wide Field 
X-ray Telescope mission (WFXT) can effec- 
tively address some of these fundamental is- 
sues through sensitive, large area X-ray sur- 
veys. 

2. The WFXT surveys at a glance 

As detailed in Rosati et al. (2010, this volume), 
WFXT is a mission thought and designed to 
perform X-ray surveys, featuring a large (1 
deg 2 ) field-of-view (FOV), large (1 m 2 at 1 
keV) effective area, and sharp (5" HEW) res- 
olution, constant across the FOV (the quoted 
values refer to the goal mission design). The 
observational survey strategy with WFXT will 
consist of three main X-ray surveys with dif- 
ferent area and depth to sample objects in a 
wide range of redshifts and luminosities. To 
illustrate the power of the WFXT surveys, a 
simple excercise can be performed by rescal- 
ing the number of X-ray sources obtained from 
well-known X-ray surveys of similar sensi- 
tivity. The 20000 deg 2 WFXT- Wide survey 
will cover about 2000 times the area of the 



XBootes survey dMurrav et al.ll2005h . in which 
more than 3000 X-ray sources have been de- 
tected, with the same depth. Similarly, the 3000 
deg 2 WFXT-Medium and 100 deg 2 WFXT- 
Deep surveys will cover 3000 times the ar ea of 
the C-COSMOS survey dElvis et al.ll2009l) and 
1000 times the a rea of the 2Ms CDFS survey 
dLuo et alJl20"08h . in which 1700 and 450 ob- 
jects have been detected, respectively. By sum- 
ming these numbers it is easy to see that the 
total source sample obtained by WFXT, mainly 
composed by AGN, will contain more than 10 
million objects. 

A more precise estimate of the number of 
AGN to be detected is obtained by considering 
the logN-logS relationships in the soft, 0.5-2 
keV, and hard, 2-7 keV, bands. About 15 mil- 
lions AGN detections are expected in the soft 
band up to z > 6 and about 4 millions in the 
hard band. Remarkably, for a very large num- 
ber of objects it would be possible to obtain 
an accurate spectral characterization over the 
0.5-7 keV WFXT band pass and derive physi- 
cal properties such as the X-ray absorption (see 
Section 4 and Matt & Bianchi, 2010, this vol- 
ume). Indeed, it is worth noting that, while be- 
ing most effective at 1 keV, the WFXT collect- 
ing area at > 4 keV is still significantly large: 
in the goal design, the effective area at 5 keV is 
equal to that of XMM (pn+2MOS combined) 
at the same energy. 

In the following Sections we will present 
a few unique science cases which can be ad- 
dressed only with wide area and sensitive X- 
ray surveys. 

3. Supermassive Black Holes in the 
early Universe 

Most stars formed at 0.5 < z < 3, when 
SMBHs were also growing most of their mass, 
but the assembly of the first organized struc- 
tures started at earlier epochs, as soon as 
baryons were able to cool within dark matter 
halos. Discovering the first galaxies and black 
holes ever formed and understanding how their 
(concurrent?) formation takes place is a fasci- 
nating subject to which more and more efforts 
will be devoted in the next decade. To date, 
several tens of galaxies and about 20 QSOs 
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have been confirme d spectroscopicaliy at red- 
shifts above 6 (e.g . lTaniguchil l2008: Fan l2006t 
IWillott et al.ll2009l). Besides a GRB discovered 
atz ~ 8.2 dSalvaterra et1ifll2009UTanvir et al. 1 
2009) and a galaxy at z = 8.6 dLehnert et al.1 



2010), no other object is known spectroscopi- 



cally at z > 7, although more than 100 galaxy 
candidates have been recently isola ted through 
deep n ear-IR imaging with WFC3 dOesch et al.l 
120101: IWilkins et all 1201 Oh . The most distant 
SMBHs known to date are three QSOs at 



6.4 discovered through the SDSS (Fan 



2006) and the CFHQS survey (Will ottet al 



2009^ 2010). These two surveys have been the 



main resource for investigating the most dis- 
tant QSOs in the past decade. About 20 QSOs 
at z > 5.7 have been detected in the SDSS main 
survey, which represent the brightest tail of the 
early QSO population (logL x - 45.5; logLboi ~ 
47). Less luminous objects have been detected 
in the SDSS deep stripe dJiang et all 120091) 
and CFHQS. The BH masses measured for 
the SDS S bright objects are of the order of 
10 9 M dKurk et al.ll2007l l2009t) . which must 
have been built in less than 1 Gyr, i.e. the 
age of the Universe at z = 6. These giant 
black holes are thought to have formed through 
mass accretion onto smaller seeds with mass 
ranging from 1O 2 M , as propo sed for the rem- 
nants of massive, PopIII stars ( Ma dau & Reesl 
1200 lb . to 1O 4 M , as proposed for the products 
of direct collapse of large molecular clouds 
dVolonteri et alJl2008l) . Whatever the seed ori- 
gin is, to account for the observed BH mass 
distribution at z ~ 6, BH growth must have pro- 
ceeded almost continuously at (or even above) 
the Eddington limit for a time interval of 
about 1 Gyr. Recent hydrodynamical simula- 
tions have shown that, within large dark mat- 
ter halos, merging between proto-galaxies at 
Z ~ 14 with seed BH masses of 10 4 M o may 
trigger Eddington limited nuclear acti vity to 
produ ce a 10 9 solar mass BH by z ~ 6 dLi et al.l 
2007). The overall frequency and efficiency by 
which this merging and fueling mechanism can 
work is however unknown. 

Because of the paucity of observational 
constraints, the study of BH and galaxy forma- 
tion at early epochs remained speculative. For 
example, the relation between the BH mass and 



the galaxy gas mass has been measured for a 
few bright QSOs at z ~ 5 -6, in which the ratio 
between the BH mass and the dynamical mass 
within a few kpc radius, as measured from CO 
line observations, is of the order of 0.02-0.1 
dWalter et al.ll2004l) . This is more than one or- 
der of magnitude larger than the ratio between 
the BH mass and stellar bulge mass measured 
in local galaxies. Despite the uncertainties on 
these measurements (see e.g. Narayana n et al.l 
120081) . the possibility that SMBHs are leading 
the formation of the bulge in proto-galaxies, 
poses severe constraints to galaxy formation 
models. Large statistical samples of z > 6 AGN 
are needed to understand the relation between 
the BH growth and formation of stars in galax- 
ies at their birth. B ase d on SAMs (see e .g. 
Hopki ns et al . 2008 and Marulli et al. 2008 for 
recent works), the luminosity function and spa- 
tial distribution (clustering) of AGN at z > 6 
would constrain: i) where, i.e. in which dark 
matter halos, they form; ii) what their average 
accretion rate is; iii) what the triggering mech- 
anism (e.g. galaxy interactions vs secular pro- 
cesses) is. 

The accreting SMBHs at z > 6 known to 
date have been all selected as /-band dropouts 
(i — z 2) in optical surveys. The sky den- 
sity of the bright {zab < 20.2) QSOs at z > 
6 in the SDSS main survey is 1/470 deg 2 . 
Fainter objects (zab ~ 21 - 22) have densi- 
ties of 1 every -30 -40 deg 2 dJiang et al.ll2009t 
IWillottetal .1120091) . The discovery of signifi- 
cant samples of objects therefore relies on the 
large areas covered by these surveys (>8000 
deg 2 for the SDSS, -900 deg 2 for the CFHQS). 
Based on the SDSS and CFHQS samples, a 
few attempts to measure the luminosity func- 
tion of z ~ 6 QSOs have been per formed, the 
most r ecen t of which are those by IJiang et all 
(2009) and lWillottetaD (l2010h . While the un- 
certainties at low luminosities (M1450 > -25) 
are still substantial, the space density of lu- 
minous QSOs (M1450 < -26.5; L boi > 10 47 
erg s" 1 ) is relatively well constrained, decreas- 
ing exponentially from z ~ 3 to z ~ 6 (see 
Fig. [TJ. Furthermore, multiwavelength studies 
of SDSS QSOs at z ~ 6 showed that these ob- 
jects have metallicities and spectral energy dis- 
tributions pretty similar to those of lower red- 
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Fig. 1. The space density of unabsorbed AGN with 
logL x = 43.5 and \ozL x = 4 5.5. Big open circles are 
from lHasinger et al.l l l2005h . Small black symbols 
have been obtained from the space density of op - 
tically luminous SDSS QSOs dRichards et al]|2006h 
by assuming a ox = -1.5. The space density of lu- 
minous QSOs is found to decrease from z ~ 3 to 
z ~ 6 and is well fitte d by an exponen tial decline 
curve (dotted line from lGilli et alj|2007h . The space 
density of lower luminosity objects is unconstrained 
for z > 4, and we show two alternative working hy- 
potheses: either a "decline" scenario, as observed 
for luminous QSOs, or a constant density model, 
which is referred to as "maximum XLF". 



shifts QSOs, suggesting their are already ma- 
ture objects in a young Universe (see however 
Jiang et afl l2010l for two counter-examples). 

Most of the multiwavelength studies men- 
tioned above refer to optical QSOs above the 
knee of the luminosity function. These likely 
represent a minor fraction of the active BH 
population at z > 6, which is expected to be 
made primarily by less massive, « 10 6 M Q ,and 
less luminous, 10 44 erg s _1 , objects. 

When compared with optically selected 
objects at the same redshift, X-ray selected 
AGN are on average less luminous, either in- 
trinsically or because of obscuration effects. 
Therefore X-ray selection might be the key to 
sample the bulk of the high-z AGN popula- 
tion. Several arguments indeed suggest that ob- 
scured AGN should be abundant even at very 



high redshift. First, observations and model- 
ing showed that obscured AGN outnumber un- 
obscured ones by a significant factor up to 
z < 4 and it is then reasonable to extrapo- 
late that even at z > 6. Second, most current 
models of galaxy formation postulate that the 
early phases of acc retion onto se ed black holes 
are obscured (e.g. Menci et aQ l2008f) . Third, 
optical/near-IR spectroscopy and IR/sub-mm 
imaging of z ~ 6 QSOs showed that dust 
and metals are abu n dant in their inner r egions 
dJuarez et alj 120091; iBeelen et alj 120061) . Dust 
and metals must have then formed in large 
quantities at z > 6 and can effectively absorb 
the nuclear radiation from the optical regime 
to the soft X-rays. 

Among the objects discovered so far at 
z > 6 no one is obscured, since these have 
been selected only through optical color selec- 
tion criteria. The most distant object discov- 
ered to date through X-ray selection is an AGN 
at z = 5.4 in the COSMOS field (Civano et al. 
in preparation) and only a few are known at 
z > 4. The lack of sufficient sky area covered 
to deep sensitivities is the reason behind that. 
The discovery space for early obscured AGN is 
therefore huge, and can have a significant im- 
pact on our understanding of BH and galaxy 
formation. 

Different routes can be tried to estimate the 
space and surface density of X-ray emitting 
high-z AGN and hence make forecasts for the 
samples to be observed in the WFXT surveys. 
As a first step, we considered extrapolations 
based on our current knowledge of the AGN 
X-ray luminosity function (XLF) and evolu- 
tion. As shown in Fig. Q] the space density 
of luminous, unobscured QSOs with M1450 < 
-26.5 has been determine d to steadily decreas e 
from z ~ 3 to z ~ 6 dRichards et alJl2006h . 
When their optical luminosity is converted to 
the X-rays (using a ox = —1.5), there is a good 
match between their space density and that 
determined by Hasinger et al.l (120051) for soft 
X-ray selected QSOs in the overlapping red- 
shift range. The behaviour of lower luminos- 
ity sources, log L A = 43.5, is instead known 
only u p to z ~ 4 dYencho et alj2009tlArr"d et alj 
|2010|) and extrapolations to higher redshifts are 
highly uncertain (see Brusa et al. 2010, this 
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Fig. 2. Left: Expected number counts in the 0.5-2 keV band for z > 6 AGN according to different models 
as labeled (see text for details). The datapoints and shaded regions show the observed, total logN-logS of 
X-ray sources. The sensitivity limits of the three WFXT surveys are shown as vertical dashed lines. Right: 
As in the left panel but for z > 8 AGN. 



volume). Starting fr om the X-ray b ackground 
synthesis model bv lGilli et al.l (120071) . we con- 
sidered a first scenario in which the space den- 
sity of AGN with logL, < 44 undergoes the 
same exponential decline as observed for lumi- 
nous QSOs - which we will call the "decline" 
scenario - and a second, perhaps extreme, sce- 
nario in which it stays constant above a red- 
shift of ~ 4, which will be referred to as 
the "maximum XLF" scenario (see Fig. 1). 
The 0.5-2 keV logN-logS of z > 6 AGN ex- 
pected in the two above mentioned scenarios is 
shown in Fig. [2] where we also show a range 
of predictions based o n different SAMs of 
BH/galaxy coevolu tion ( Salvat erra et a J 20071: 



2008). 



iRhook & Haehnellll2008t iMarulli et al 
The limiting fluxes of the three WFXT sur- 
veys are also shown. It is clear how the z > 6 
AGN Universe is a completely uncharted terri- 
tory, the various predictions differing by a few 
orders of magnitude already at fluxes around 
10' 16 cgs. The predic tion s by the SAMs of 



ISalvaterra etail d2007l) and lRhook & Haehneltl 
(2008) are even more optimistic than the "max- 
imum XLF" scenario, predicting about 500- 
600 z > 6 AGN deg 2 at the WFXT-De e p lim- 
iting flux. The model bv lMarulli etafl 00081 



is instead the most pessimitic, with 15 z > 6 
AGN deg~ 2 at the same limiting flux. The pre- 
dictions by SAMs depend on a pretty large 
number of parameters such as the mass of seed 
BHs, their location in the density field (i.e. 
the abundance of the dark matter halos host- 
ing them), the recipes used for accretion, the 
AGN lightcurve. Although different assump- 
tions are made by different models, the ~ 3 
orders-of-magnitude difference at faint X-ray 



fluxes between the predictions by Marulli et al. 



(2008) on the one hand, and Salvaterr a et al. 



(120071) and IRhook & Haehneltl d2008l) on the 
other hand, seems to be primarily related to 
the ass umed space d e nsity o f seed black holes. 
In the IMarulli et all d2008l) model, seed BHs 
are placed in each halo that can be resolved 
by the Millennium simulation, i.e. in those 
ha los with mass above ~ 10 10 M , while 
in ISalvaterra et al.l d2007) BHs populate mini- 
halos with mass as low as ~ 10 7 ~ 8 M Q , 
which are therefore much more abundant. In 
IRhook & Haehneltl (12008) BHs are assumed to 
radiate at the Eddington limit and their mass 
to be proportional to that of the hosting halos. 
This implies that, if small mass seed are as- 
sumed, these populate low mass, abundant ha- 
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Fig. 3. Left: Number of AGN at z > 6 expected from the "decline" model (see text) for different combina- 
tion of survey area vs 0.5-2 keV limiting flux. Dotted lines are the locii of equal AGN number as labeled 
(labels are in log units). Small black squares show the major Chandra and XMM surveys (either performed 
or ongoing, including the 4Ms CDFS, the 2Ms CDFN, C-COSMOS, AEGIS, XBootes and many others). 
Big black squares show the current coverage of the Chandra (CSC) and XMM (2XMM-pn) archives. The 
deep and wide surveys to be performed with eROSITA in 4 years are shown as magenta circles (based 
on Cappelluti et al. 2010, this volume). The three WFXT surveys (5 year program) are shown as rotated 
squares. Less than one z > 6 AGN is expected in each of the current X-ray surveys. eROSITA will return 
about 40 objects with the current design. WFXT will detect about 1600 objects. 



los. If large mass seeds are assumed, these pop- 
ulate less abundant halos but their luminosity, 
and hence detectability, is higher. In each case, 
surface densities as high as 500 deg~ 2 at the 
WFXT-Deep limiting flux are reached. 

Observations of significant samples at z > 
6 would constrain the physics of early BH for- 
mation disentangling between these scenarios. 
The surveys performed by WFXT have the 
power to provide such large statistical samples. 
As a reference model, we will consider the "de- 
cline" model outlined above, which is found to 
be in excellent agreement with th e 0.5-2 keV 
logN-logS of AGN at z > 3,4,5 (iBrusa et all 



I2009L Civano et al. in prep.). About 1600 AGN 
at z > 6 are expected to be detected by WFXT, 
and about 70 ot them should be at z > 8. A 
summary of the predictions for AGN at z > 6 
and z > 8 in the three WFXT surveys for the 
decline and maximum XLF scenarios is shown 
in Table 1 . In the WFXT-Deep survey it will be 
possible to detect in significant numbers z > 6 
AGN down to logL v = 43.1 and z > 8 AGN 
down to logLj = 43.4, making possible to de- 
termine the shape of the high-z XLF. To com- 
pare the high-z AGN yields from WFXT sur- 
veys with what is expected from major cur- 
rent X-ray surveys and with other X-ray mis- 
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Table 1. Number and minimum 0.5-2 keV 
luminosity of high-z AGN expected in the 
WFXT surveys according to the two evolution 
models ("decline" and "max LF") described in 
the text. 



Quantity 




Survey 






Deep 


Medium 


Wide 


z > 6 


logL"""(0.5-2 keV) 


43.1 


44.3 


45.1 


N. AGN (decline) 


300 


1000 


300 


N. AGN (max LF) 


15000 


2300 


300 


z> 8 


logL"""(0.5-2 keV) 


43.4 


44.6 


45.4 


N. AGN (decline) 


20 


45 


10 


N. AGN (max LF) 


4300 


210 


10 



sions which are either proposed or planned, we 
show in Fig. [3] the number of z > 6 AGN to 
be detected with different combination of sur- 
veyed sky area and soft limiting flux according 
to the "decline" scenario. It is evident that even 
the major current X-ray surveys (e.g. CDFS, 
CDFN, COSMOS, AEGIS, X-Bootes) should 
return less than one object per field. In the full 
Chandra and XMM archives (labeled as CSC 
and 2XMM-pn, respectively) one would expect 
less than 10 objects in total. With the current 
mission design and survey strategy (Cappelluti 
et al. 2010, this volume), the eROSITA satel- 
lite, planned for launch in 2012, is expected 
to return about 40 z > 6 AGN, most of them 
at the luminosities of bright SDSS QSOs. The 
big leap is clearly expected to be performed 
by WFXT. The mission design and observ- 
ing strategy studied for the Internation X-ray 
Observatory (IXO) are still uncertain: IXO will 
probably return a sample of z > 6 AGN with 
size in between those by eROSITA and WFXT 
(see Comastri et al. 2010, this volume). 

The identification of these objects will 
clearly require wide area surveys with deep 
multiband optical and near-IR imaging like 
e.g. the LSST surveys (Brusa et al. 2010, this 
volume). On the other hand, the WFXT sur- 
veys will represent a perfect complement for 
all optical and near-IR campaigns that search 
the unobscured part of the high-z AGN pop- 
ulation as image dropouts. Late type brown 
dwarfs are the main contaminants in optically 



selected samples, being 15 times more abun- 
dant t han high-z QSO s of comparable magni- 
tudes dFan et alj|200ll) . Even when using near- 
IR colors to separate them from high-z QSO 
candidates, the success ra te of optical spec - 
troscopy is only 20-30% dJiang et alj l2009h . 
If sensitive X-ray images were available un- 
derneath each optical dropout, a detection in 
the X-rays would almost automatically ensure 
that the object is an AGN, since brown dwarfs 
of comparable optical magnitudes are ~ 300 
times fainter in the X-rays. LSST and WFXT 
surveys will then reinforce each other in the 
search of z > 6 AGN. 



4. Evolution of the obscuration 

A built-in feature of the BH/galaxy evolution- 
ary sequence described in the Introduction is 
that an obscured accretion phase preceeds a 
clean accretion phase, at least in powerful, 
QSO-like objects. One may therefore wonder 
whether the fraction of obscured AGN was 
higher in the past. This indeed depends on 
many parameters such as i) the physical scale 
of the absorbing gas and how this is driven to- 
wards the BH; ii) the relative timescales of the 
obscured and unobscured phases; iii) whether 
the absorbed-to-unabsorbed AGN transition 
occurs also in low-mass/low-luminosity ob- 
jects (i.e. Seyfert galaxies). Not many theo- 
retical predictions on the evolution of the ob- 
scured AGN fraction are available in the lit- 
erature. Models that relate the obscuration on 
nuclear scales to the availability of gas in the 
host galaxy generally predict an increase of 
the obscured AG N fraction with redshift (e.g 
Menci et al. 2008) since the gas mass in galax- 
ies was larger in the past. Some others mod- 
els, which anti-correlate the covering factor of 
the obscuring medium to the BH mass and 
then follow the evolution of the BH mass func- 
tion u sing empirical relations (Lamastr a et al] 
2008), however do not predict such an in- 
crease. The situation is also debated from 
an observational point of view. An increase 
of the obscured AGN fraction with redshift 
am ong X-ray selecte d AGN h as been observed 
by lLa Franca et all (120051): iTreister & Urrvl 
(2006); Hasingerl d2008l) : iTrump et al l (12009). 
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Fig. 4. Left: Redshift distribution of the 300,000 AGN with good X-ray spectra (i.e. with more than 1000 
photons) in the three WFXT surveys. The Deep survey is providing most of the good X-ray spectra at z > 1. 
Right: As in the left panel but splitting the total sample of 300,000 AGN into unabsorbed (logN H < 22) and 
absorbed (logNn > 22) objects. 



but other works did not find any evidence of 
this tr end dUeda et alJ 120031: iDwellv & Page] 
120061: iGilli et alJ 120071) . Selection effects in 
these computations are very important and 
could mimic a spurious evoluti on of the 
obscur ed AGN fraction (see e.g. IGilli et ail 
I2010ah . 

Much of this uncertainty can be related 
to the lack of large statistical samples of 
AGN with high-quality X-ray spectra. Using 
samples of limited size, which cannot be 
split into narrow luminosity and redshift bins, 
it is often complicated to disentangle any 
redshift-dependent from luminosity-dependent 
effects (it is generally agreed that the ob- 
scured AGN fraction decreases with lumi- 
nosity, but the slope and normalization of 
this relation differ from paper to paper). 
Furthermore, the absorption column density is 
often estimated either from the X-ray hard- 
ness ratio, or from the absence of broad 
emission lines in the optical spectrum. Both 
methods suffer fro m caveats and limitations 
(iBrusa et alJ 120101) . and ideally the absorb- 
ing column should be measured through X- 
ray spectra with good photon statistics. The 
most recent XLF determinations contain < 



2000 AGN dHasingerl 120081: ISilverman et 
2008: lYencho et alj|2009t lEbrero et al.ll200' 



al. 

)9; 



Aird et all 12010). most of which are just 



slightly more than simple X-ray detections. 
With WFXT one would expect to detect about 
300,000 AGN with good X-ray spectra, i.e. 
with more than 1000 photons each, in a broad 
redshift range (z ~ - 5). As shown in 
Fig.|4]/e/f, the WFXT-Deep survey is providing 
most of the good X-ray spectra at z > 1 , while 
the Medium and Wide surveys provide the 
largest contribution at lower redshifts. About 
1/3 of the detected objects with good X-ray 
spectra are expected to be absorbed by column 
densities above 10 22 cm~ 2 (see Figj4] right). 
Future wide area optical and near-IR surveys 
should provide redshifts, either spectroscopic 
or photometric, for a significant fraction of 
these sources, making possible to divide the 
sample into fine redshift, luminosity and ob- 
scuration bins and therefore map the cosmo- 
logical evolution of nuclear obscuration. 



5. The most obscured AGN 

If the overall evolution of obscured AGN is 
uncertain, the evolution of the most heav- 
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2 ^ 5 2 1 2 5 

channel energy (keV) channel energy (keV) 

Fig. 5. Simulated WFXT spectra of CDFS-202 and CDFS-153, two high-z CT AGN with / 2 _ 10 > 10" 15 
erg cm -2 s _I observed in the XMM-CDFS: the goal effective area and an exposure time of 400 ks were 
assumed. No background is included, w hich however is not e xpected to alter significantly the S/N ratio (see 
text).The 100 deg 2 WFXT-Deep survey dMurrav et alJl2008L Rosati et al. 2010, this volume) is expected to 
reveal about 500 objects like these at z > 1 in terms of obscuration and photon statistics. 



ily obscured and elusive ones, the so-called 
Compton-Thick AGN (Nh > 10 24 ctrT 2 ; here- 
after CT AGN), is completely unknown, mak- 
ing the census of accreting black holes largely 
incomplete. 

Only -50 bona-fi.de CT AGN, i.e. certified 
by X-ray spe ctral analysis, are known in the lo- 
cal U niverse ((Comastri 2 0041: iDella Ceca et al.1 

2008) , but their abundance has nonetheless 
been estimated to be comparable to that of less 
obscured ones (Risa litietali ri999). At higher 
redshifts, an integrated constrain to the space 
density of CT AGN can be obtained from the 
spectrum of the X-ray background. Depending 
on the assumptions, synthesis models of the 
XRB predict that from -10 to -30% of the 
XRB peak emission at 30 keV is produced 
by the emissio n of CT AGN integrated over 
all re dshifts dGilli et al.1 120071: iTreister et afl 

2009) , but the main contribution likely arises 
at z - 1 • In the absence of any information, the 
luminosity function and evolution of CT AGN 
in XRB synthesis models have been usually as- 
sumed to be equal to those of less obscured 
ones. Because of the uncertainties in the aver- 
age spectrum of CT AGN and in their evolu- 
tion, however, the constraints to the CT space 
density from the XRB spectrum remain rather 
loose. 



In recent years there have been many at- 
tempts to constrain the space density of CT 
AGN in different luminosity and redshift inter- 
vals exploiting different selection techniques. 
Very hard (> 10 keV) X-ray surveys are still 
limited in sensit ivity and are just sa mpling the 
local Universe (Ma lizia et alj|2009l) . To sam- 
ple AGN at higher redshifts, and in particu- 
lar at z - 1, one needs to rely on deep X-ray 
surveys in the 2-10 keV band and try to se- 
lect CT AGN either directly from X-ray spec- 
troscopy, or by comparing the measured, ob- 
scured X-ray emission (if any) with some other 
indicator of the intrinsic nuclear power: IR se- 
lection can track the nuclear emission as re- 
processed by the dusty absorber; narrow opti- 
cal emission lines can sample the gas ionized 
by the nucleus on scales free from obscuration. 
X-ray stacking of IR-selected sources not indi- 
vidually detected in the X-rays has been used 
to estimate the space density of CT AGN at 
z - 1-2 dDaddi et al Il2007t lAlexander et all 
l2008tlFiore et al.l2009tlBauer et al.l2010h . The 
comparison between the [O III]5007 and X-ray 
flux has been used to select X-ray underlumi- 
nous QSOs and then estimate the densit y of 
CT QSOs at z = 0.5 dVignali et al.ll2010h . To 
sample the popula tion of CT around z - 1 , 
iGilli et all (1201 Obi) recently devised a selec- 
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tion method based on the [Ne V]3426 emission 
line, which can be applied to optical spectro- 
scopic surveys with deep X-ray coverage and 
seems to deliver clean, albeit not complete, 
samples of CT objects at z > 0.8, i.e. at red- 
shifts not reachable with [O III]5007 selection. 
Despite these numerous efforts, [Ne V], [O III], 
and IR selection are all indirect ways to select 
CT AGN, since the CT nature of an object is 
simply inferred from the faintness of its X-ray 
emission relative to an indicator of the intrin- 
sic power, and therefore may suffer from severe 
systematic uncertainties. 

To unambiguously obtain samples of bona- 
fi.de CT AGN over a broad redshift range, very 
deep X-ray exposures are needed such as the 
2.5 Ms XMM survey in the CDFS (Comastri 
et al. 2010, in prep.). A number of CT can- 
didates, identified in the 1 Ms CDFS catalog 
on th e basis of their fl at (low quality) spec- 
trum dTozzi et al.ll2006h . are indeed being con- 
firmed as such by the higher quality XMM 
spectra (Comastri et al. 2010, in prep.), includ- 
ing the well-known CT candi date CDFS -202 at 
z = 3.7 dNorman et al.ll2002l) . Because of their 
limited sky coverage, however, only a few tens 
of bona-fide high-z CT AGN are expected to be 
detected in current deep X-ray surveys, making 
population studies of CT AGN problematic. 

The WFXT surveys are expected to deter- 
mine the cosmological evolution of bona-fide 
CT AGN up to z ~ 3. Based o n the synthe- 
sis model by iGilli et alj d2007l) . the WFXT- 
Deep survey will return a sample of ~ 500 ob- 
jects at z > 1 which are bona-fide CT AGN, 
i.e. with more than 500 net counts in the 0.5- 
7 keV band (the number of simple detections 
of CT objects will be obviously much larger). 
In Fig. |5]we show two such objects (CDF-202 
and CDFS- 153, another CT AGN, at z = 1.53, 
found in the XMM-CDFS) simulated using the 
WFXT response matrices for the goal design. 
No background is assumed in the simulation, 
but, based on the level estimated for low-earth 
orbits (see Ettori & Molendi 2010, this vol- 
ume), only about 20 background photons are 
expected above 1 keV, which would therefore 
not alter significantly the spectral quality. It is 
evident that 500 X-ray photons are sufficient to 
unambiguously reveal their CT nature: in prin- 



ciple, based on the iron Ka line, it would be 
possible to determine also their redshift with- 
out the need of optical spectroscopy. In to- 
tal, 500, 270, 60 and 12 bona-fide CT AGN 
at redshifts above 1, 2, 3 and 4, respectively, 
are expected in the WFXT-Deep survey. Future 
missions sensitive to energies above 10 keV, 
such as NuSTAR and ASTRO-H (approved by 
NASA and JAXA, respectively) or EXIST and 
NHXM (proposed to NASA and ASI, respec- 
tively) are expected to allow population studies 
of CT AGN at z < 1, with a peak in the red- 
shift distribution at z ~ 0.3 - 0.4. The WFXT 
mission appears to uniquely complement these 
and extend them to z > 1. Only IXO, on a 
longer timescale and depending on the adopted 
survey strategy and mission design, could pos- 
sibly provide samples of distant CT objects 
matching in size those expected from WFXT. 



6. Conclusions 

About 15 millions of AGN in a very broad 
redshift and luminosity range will be detected 
by the WFXT surveys. Some of the major 
issues related to the evolution of accreting 
SMBHs are expected to be solved by this 
sample as detailed below. 

• WFXT will break through the high-z 
Universe: more than 1600 AGN will be 
observed at z > 6, allowing population studies 
at these redshifts and providing an unvaluable 
complement for future wide area optical and 
near-IR surveys searching for black holes at 
the highest redshifts. Such a large sample is 
a unique feature of the WFXT surveys and 
cannot be matched by any other planned or 
proposed X-ray mission. 

• Good X-ray spectra, with more than 1000 
photons, will be obtained for 300,000 AGN 
allowing accurate measurements of their 
absorbing column density. This will make 
possible for the first time to measure the 
evolution of nuclear obscuration with cosmic 
time up to z ~ 5 and verify its connection with 
e.g. star formation. 
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• About 500 bona fide, i.e. certified by X- 
ray spectral analysis, CT AGN at z > 1 will 
be found in the WFXT-Deep survey. This will 
complement population studies at lower red- 
shifts obtained by future high-energy surveys 
such as those performed by the approved mis- 
sions NuSTAR and Astro-H and will make 
possible to detemine the abundance and evo- 
lution of this still missing BH population. 
Performing population studies of distant heav- 
ily obscured objects, and determining their 
relevance in the census of accreting black 
holes and evolutionaly path of galaxies is an- 
other unique science case to be carried out by 
WFXT. 
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